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ABSTRACT 

We use a set of magnetohydrodynamics (MHD) simulations of fully-developed (driven) turbulence 
to study the anisotropy in the velocity field that is induced by the presence of the magnetic field. 
In our models we study turbulence characterized by sonic Mach numbers Ms from 0.7 to 7.5, and 
Alfven Mach numbers from 0.4 to 7.7. These are used to produce synthetic observations (centroid 
maps) that are analyzed. To study the effect of large scale density fluctuations and of white noise 
we have modified the density fields and obtained new centroid maps, which are analyzed. We show 
that restricting the range of scales at which the anisotropy is measured makes the method robust 
against such fluctuations. We show that the anisotropy in the structure function of the maps reveals 
the direction of the magnetic field for Ma ^ 1-5, regardless of the sonic Mach number. We found 
that the degree of anisotropy can be used to determine the degree of magnetization (i.e. Ma) for 
Ma ^ 1.5. To do this, one needs an additional measure of the sonic Mach number and an estimate 
of the LOS magnetic field, both feasible by other techniques, offering a new opportunity to study the 
magnetization state of the interstellar medium. 

Subject headings: ISM: general — ISM: structure — magnetohydrodynamics (MHD) — radio lines: 
ISM — turbulence 



1. INTRODUCTION 

It is well accepted that turbulence plays a central 
role in the dynamics and transport phenomena in the 
interstellar medium (ISM). For an overview we refer 
the re a der to the recent reviews by lElmcgrc en fc Scalo 
(20041); pac Low fc Klessen' (2004) ; Ballesteros-Paredes 



et al. 



(2DD7); McKee & Ostriker (2007), and references 



therein. Moreover, turbulence in the iSM is magnetized, 
so in order to understand its properties one has to ad- 
dress those of the magnetic field as well. 

The presence of a magnetic field introduces a prefer- 
ential direction of motion for the charged particles and 
in consequence makes the turbulent cascade anisotropic, 
which has been known for some time now ( [Montgomery! 
[1982; Higdon 1984). In a turbulent magnetized medium 
the kinetic energy of large-scale motions is larger than 
those at small-scales, but the local magnetic field (thus 
the magnetic energy) is comparable at all scales. There- 
fore at large-scales the magnetic field does not dom- 
inate dynamically, but it becomes more important as 
we go into smaller scales. The result are elongated ed- 
dies, which become more elonga ted as the energy cas- 
cade goes down to smaller scales (|Cho fc Vishniac||2000 



Maron fc Goldreich 
2002' henceforth ICL 



2001 |Cho, Lazarian, mac 
'02p. Recent decades have seen im- 



portant progress in the understanding of the magnetic 
turbulence. A pioneering study is the self-consistent 
model of magnetohy d rodyn amic (M HD) tu rbulence of 
Goldreich fc Sridhar[ ([19951 hereafter [GS95[ ) where the 



concept of a scale dependent anisotropy is built ir{^ The 
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'I'he initial publication did not have tHe^concept or the local 
magnetic field built in. In fact, the closure relations in the paper 



|GS95| model w as later supported by the numerical sim- 
ulation s of Ch o fc Vishniacj ( [2000| ; [Maron fc Goldreich 
(2001 1; |CLV()2) . I'he ideas on compressible turbulence 
sketched in |GS95[ have proved profound, and have been 
supported and developed in further studies ([L ithwick & 
Goldreichl200Tl|Cho fc Lazarian |20031 IKowafe 



owal"fc Lazarian 



2010[ ). while some aspects of GS95 scaling were claimed 
to be controversial and not consistent with numerical 
simulations, e.g. the predicted Kol mogorov-type index 



ian 


2010 


cal 


simui 



have revealed the limitations ot the numeri- 
ations which induced these controversies. All 



in all, we believe that the GS95 model (with the adjust- 
ments and improvements introduced in subsequent pub- 
lications) provides the best representation of the MHD 
turbulence statistics. 



The GS95 model is w ell described in many reviews (see 
Lazarian fc Cho|[2005 ) and we do not dwell upon its de- 



tails. We just mention that the qualitative understand- 
ing of the model can be obtained if one assumes that 
mixing motions perpendicular to the local magnetic field 
create a turbulent Kolmogorov-like cascade. The eddies 
are elongated along the local direction of magnetic field 
and the relation between the parallel and perpendicular 
scales of the eddies is given by the so-called "critically 
balance," which is refiected in the equality of the eddy 
turnover time and the timescale of propagation of Alfven 
waves along the magnetic field of the eddy. The tensor 
describing the turbulent magnetic field is presented in 



were written in the mean field frame of reference. The correct 
understanding that the eddies can be described only in the system 
of reference related to the lo cal magnetic field of the eddies in 
question was obt ained later (Lazarian & Vishniac||1999| ICho &| 
[Vishniac_2000 , Maron fc Goldrcich^2001 . ULV02| ) 
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The underlying anisotropies of turbulence result in 
the anisotropies of the observed statistics of turbulence. 
However, the anisotropy of interstellar turbulence which 
is accessible to the observer averaging emission along the 
line of sight crossing a turbulent volume is different from 
the |GS 95 predictions. In the "global" system of refer- 
ence, related to the mean magnetic field, the anisotropy 
is no longer-scale dependent but is determined by the 
anisotropy of the largest eddies. A suggestion of using 
such anisotropy to study turbulence and the direction of 
the mean rnagnet ic field was made [ Lazarian, Pogosyan^ 



fc Esquivel] ( |2002| henceforth |LPE02p , where the feasi- 
bility ot such stuaies was illustrated with synthetic spec- 
tral line emission maps obtained via MHD turbulence 
simulations. The anisotropy is readily evident from two- 
point statistics (e.g. correlation/structure functions and 
power-spectra), where more power (larger dispersion) is 
concentrated in the directi on perpendicular to the mag- 
netic field ( see for instance Esquivel et al.|[2003 Vestuto 
3J). ™ " ■ 



et al. 2003|. The aforementioned studies provided the 
framew ork for the observational studies in |Heyer et al. 



( |2008 ). 

in this paper we revisit the velocity anisotropy from nu- 
merical simulations and synthetic observations. In par- 
ticular, we study how such anisotropy depends on the 
global Alfvenic and sonic Mach numbers of the turbu- 
lence. 

If one knows the Alfven Mach number, then by know- 
ing the turbulent velocity dispersion, an estimation of 
the media magnetization (i.e. its Alfven speed) is pos- 
sible. The determination of the Alfen speed, which is a 
crucial parameter for theoretical and numerical models, 
has been a great challenge for current techniques that 
probe the ISM. 

We describe our MHD simulations in §2, and the way 
of measuring velocity anisotropy in observations in §3. 
Results of our study, namely the degree of anisotropy 
obtained as a function of Alfvenic and sonic Mach num- 
bers are presented in §4. We finish with a discussion and 
provide a summary in §5. 

2. MHD MODELS 

We use a set of three-dimensional MHD simulations of 
fully-developed (driven) turbulence to produce maps of 
velocity centroids and of the average LOS velocitjQ 

The simulations presented in this work were obtained 
solving the ideal MHD equations in a periodic box. 



dpv 

~dt 



pvv 



~dt 



dp 
dt 

I 



V- (pv) ^ 0, 



An 



BB 



(1) 
(2) 



Vx(vxB) = 0, (3) 



with an isothermal equation of state {p — c^p, where p is 
the gas pressure, Cg the sound speed, and p the mass den- 

^ The average LOS velocity traces the statistics of velocity in the 
same manner that column density traces density. However, it can 
not be obtained directly from observations, while column density 
can. In previous works we have ter med it "inte grated velocity" 
(e.g. [Lazarian fc Esquivcl_2003 Esquivel fc Lazarian^2005j l. 



TABLE 1 

Parameters of the MHD simulations. 



Model 


VAX) 


( Pn ,^ « n ) 
\ y " ^ 1 ^' / 


Ms 


Ma 


Resolution 


Ml 


0.1 


0.01 


~ 7.6 


~ 7.6 


512=* 


M2 


0.1 


0.1 


~ 2.4 


~ 7.7 


512=* 


M3 


0.1 


1.0 


~ 0.8 


~ 7.7 


5123 


M4 


0.5 


0.01 


~ 7.5 


~ 1.5 


256^ 


M5 


0.5 


0.1 


~ 2.3 


~ 1.5 


256^ 


M6 


0.5 


1.0 


~ 0.7 


~ 1.4 


256=* 


M7 


1.0 


0.01 


~ 7.4 


~ 0.7 


512^ 


M8 


1.0 


0.1 


~ 2.3 


~ 0.7 


512^ 


M9 


1.0 


1.0 


~ 0.7 


~ 0.7 


512^ 


MIO 


2.5 


0.01 


~ 9.8 


~ 0.4 


256^ 


Mil 


2.5 


0.1 


~ 3.3 


~ 0.4 


256=* 


M12 


2.5 


1.0 


~ 1.1 


~ 0.4 


256=* 



sity) and the additional constraint of V • B = 0, achieved 
with a con strained transport (CT) algorithm (see e.g. 
Toth 2000). The integration method is a second-order- 
accu rate hybrid essentially n onoscillatory (ENO) scheme 
(see Cho & Lazarian 2002). In order to avoid spuri- 
ous ossc ilations, the code sw itches from a ENO weighted 
scheme ( Jiang fc Wu||199^ where variables are smooth) 
to a convex EiNO scheme (Shu fc Osher 1989 where 
strong gradients are found). 'I'he time is marched with 
a two stage Runge-Kutta method. The source term 
at the right hand side of equation ([2]) is a large-scale 
pseudo-random driving. Such driving is purely solenoidal 
and is performed in Fourier space at a fixed wave num- 
ber k — 2.5 (a scale of 1/2.5 of the co mputational do- 
main) . For details about the dr iving see Cho fc Lazarian 
(2003|) and |Kowal et al.| ( |2007[ ). Initially p 
Alfven speed va 



1 and the 

|B|/y'47ryO depends on the model. In 
stationary state the rms velocity is also close to unity 
('^rms ~ 0.7). The parameters that define the experi- 
ments used are the sonic and Alfvenic Mach numbers at 
the injection scale = {Vl/c^) and = {Vl/va)^ 
respectively, where Vl = Wrms is the turbulent veloc- 
ity at the injection scale, and (. . .) stands for average 
over the entire computational domain. These parame- 
ters can be controlled by the values of the gas pressure 
and the Alvfen speed at the beginning of the simulation, 
{Pgasfi) and VAfij respectively. In Table [l] we list the nu- 
merical experiments that cover sub-sonic and super-sonic 
regimes, combined with different intensities of magnetic 
field (yielding sub-Alfvenic and super-Alfvenic regimes 
as well). We have to mention that the value of Vrms in 
stationary state varies from model to model, the Mach 
numbers presented in the table are measured from the 
output of the simulations. 

The initial magnetic field is of the form B = Boxt + b, 
a uniform field Boxt plus a fluctuating part b. Initially, 
b = and Bcxt is aligned in the x— direction. When the 
simulations reach a stationary state the magnitude of the 
mean and the fluctuating parts are of the same order, 
while the mean magnetic field remains aligned with the 
x— axis. 

3. VELOCITY ANISOTROPY 

We have taken the results from the 3D-MHD simula- 
tions and produced 2D maps of the mean velocity and 
centroids of velocity, integrating along each of the car- 
dinal {x ,y ,z) directions. We consider an isothermal. 



Velocity anisotropy 
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optically thin medium, with an emissivity linearly pro- 
portional to the density (e.g. cold HI). If we integrate 
along the x— axis we can obtain a 2D map of the mean 
velocity perpendicular to the y — z plane: 



1 



Vxix,y,z)dx, 



(4) 



where is the number of cells used to discretize x, 
and Vx is the projection of the velocity field along k. 
This mean LOS velocity traces the velocity structure in a 
similar m anner as the colum n density follows the density 
structure (Esquivel & Lazarian 2005, henceforth EL05). 
However, it cannot be obtained directly from observa- 
tions. 

In real observations we are faced with a density 
weighted mean (i.e. the velocity centroids). For the 
centroid maps we used their conventional (normalized) 
form, under the assumption of an optically thin medium 



whose emissivity is p roportional to the density (see |EL05 
Esquivel et aI1[2007l ): 



Cx{y,z) 



(5) 



Jvx{x,y,z) p{x,y,z)dx 
J p{x,y,z)dx 

Analogous expressions to equations Q and ([5| can be 
used to obtain the mean velocity, or velocity centroids 
with the LOS aligned with or z. 

The two-point, second-order structure function of a 
quantity /(x) is defined as: 



5i^(r)-([/(x)-/(x + r)]' 



(6) 



where (. . . 
space (x). 



denotes an ensemble average over all 
and r is the separation or "lag". A 
closely related measure is the correlation function 
CF{r) = (/(x)/(x + r)), which differs from t he SF ba- 
sically by a constant (see for instance EL05). For the 
sake of simplicity we will use the following notation: 
5'JV,x(R-), and S'i^c,x(R.) denote the structure function 
applied to a map of mean velocity (eq. El), and a map 
of velocity centroids (eq. [s]), respectively. The lag R 
is written in upper case letters to denote that is a two- 
dimensional vector, while the x in the sub-index indicates 
that the mean velocity or centroids were obtained inte- 
grating along the x— axis [in which case R = {y, z)]. The 
power-spectrum is another recurrent tool in turbulence 
studies. It is the Fourier transform of the correlation 
function, thus it provides equivalent information. 

If the turbulence were isotropic these two point statis- 
tics would (statistically) not depend on the direction 
of the lag, or on the direction of the wave number for 
the power-spectrum. In fact it is customary to assume 
isotropy and average over all angles to present just a ra- 
dial dependence of SF(r). It is well known, however, 
that the magnetic fields breaks t he isotropy an d the tur- 



spectral index of velocity (Lazarian & Esquivel 2003 
ELOS) [Esquivel et al.||200'7[ ). 

In order to further study the impact of density fluctu- 
ations in the anisotropy that can be observed we study 
four different maps obtained for each of the simulations 
of Table [T] Firstly, a map of the mean velocity obtained 
as in equation [4] which have information exclusively of 
velocity, but can not be obtained from observations. Sec- 
ondly, a map of velocity centroids (eq. [s]), obtained 
with the density field as obtained from the simulations. 
Thirdly, another map of velocity centroids, but in this 
case modulating the density by an profile. The di- 
rection of the gradient is (x, y, z), so that it is oblique to 
both the LOS and the mean B field. To avoid the singu- 
larity at r = we have placed the origin half a pixel out- 
side the computational domain, and the resulting density 
field was rescaled to have a mean value of 1. This gradi- 
ent introduce a large-scale density variation, which could 
be encountered in observations of the ISM, but which is 
a not a product of turbulence (e.g. self-gravity). And 
fourthly, another map of centroids was obtained with a 
density field to which we have added white noise. We 
model the noise using fractional Brownian motion (fBm) 
fractal structure s, that can be characterized by a powe r 
spectrum index ( Stutzki et al.|1998 Bensch et al.|200l" ), 
which for white noise we have set to zero (flat power 
spectrum). The resulting data cubes have a Gaussian 
probability distribution function. Their dispersion have 
been scaled to the same value of the original density for 
each of the models. The noise is added to the density, 
and to preserve the density positive defined we have set 
a fioor value of 0.01 (the mean density is 1.0). 

The resulting density fields are illustrated in Fig fl] 
where we present map of the x?/— midplane cut of the 
density of one of the models (M8, super-sonic and sub- 
Alfvenic, see Table [l]). The first panel (a) displays the 
original density, the middle panel (b) shows the effect of 
the large scale gradient, and in (c) the addition of white 
noise is evident. 

In the following section we will study the correlation 
and structure functions of the maps of mean velocity and 
of centroids, and how their anisotropy depends on the 
parameters of the models. 

4. RESULTS 

Observations sample the entire LOS, and at a given 
velocity one has contribution of material that could be 
anywhere along that LOS. In some sense, one can say 
that observations average the information in the position 
perpendicular to the plane of the sky. Thus, from an 
observational point of view, it is more natural to study 
the anisotropy in a global frame, namely, the anisotropy 
with respect to the direction of the mean magnetic field 
(as op posed to the local magnetic field, from which a 
GS95 | scaling is ret rieved). 



We have shown (LPE02 Esquivel et al. 2003 EL05I 



rop 

bulence beco mes anisotropic (see Higdon|1984 and refer ^ 

ences therein, |GS95[ |Cho fc Lazarian|20(J5| tor a review) . that in such a global frame' indeed the statistics of' veloc' 

ity centroids reveal the direction of the mean magnetic 
field. In Figure [2] we present contours of equal correla- 
tion in one of the models (MS, the same used in Figure 
[T]). The mean LOS contours are remarkably similar to 
those of velocity centroids. In the figure we only show the 
results for the centroids with the original density because 
they are very similar to those obtained if the cx den- 



3.1. Additional Density Fluctuations 

Spectroscopic observations are sensitive to density and 
velocity fiuctuations simultaneously. For the particu- 
lar case of centroid maps, it was clear for instance that 
the strong fluctuations in highly supersonic turbulence 
{Ms > 2.5) affected severely our ability to determine the 
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Fig. 1. — Density cuts in the Xy—midplane. (a) Original density, 
(b) density modulated by an oc gradient (decreasing from the 
left to the right, and from the bottom to the top), and (c) the 
density with additional white noise. All the plots are in the same 
(logarithmic) scale, as indicated by the bar on the right. 



sity gradient, or white noise (see ^3.1) are included. 

In panels (a) and (c) in Figure [2] the line of sight is 
in the direction parallel to the mean magnetic field (x), 
thus the two axes shown are perpendicular to it, and the 
contours are more or less circular (isotropic). Panels (b) 
and (d) have been integrated in a direction perpendicular 
to the mean plane (z) and the resulting correlations show 
a clear anisotropy in the direction of the B field, the same 
result is obtained if we integrate along the y— axis. This 
result about the direction of the mean magnetic field (or 
more accurately, its projection onto the plane of the sky) 
is quite robust. 

To address the dependence of the anisotropy with scale, 
we start by defining a simple measure of the degree of 
anisotropy: the ratio of the structure function in the 
two directions that are perpendicular to the LOS (e.g. 
SFc^z{x, 0)/ SFc^z{0, 2/)), which we have computed for all 
the models. In Figure [3] we present an example obtained 
from the same model of the previous two Figures. In the 
top three panels (a-c) we show the results of the mean 
LOS velocity, and in the bottom panels (d-f) the results 
with the different velocity centroids. 

If this ratio is one the structure function is isotropic, 
as in the case of panels (a) and (d), which correspond 
to the structure functions when the mean magnetic field 
and the LOS are aligned. The rest of the panels are 
clearly anisotropic, with a degree of anisotropy, whose ex- 
act value depended on model. In Figure [3] we indicated 
with vertical lines the scale length of injection (turbu- 
lence forcing) with a dotted line. At such large scales, 
and down to about 1/5 of the box size (marked with 
a vertical dashed line) the effect of the forcing is evi- 
dent. For the mean LOS velocity, and the centroids with 
the original density field, we found that the anisotropy 
was virtually scale independent from the small scales up 
to separations on the order 1/5 of the computational 
box (half the size of the injection scale). For centroids 
with additional density fluctuations (dashed and dotted 
lines) we see that the anisotropy shows a flat (scale in- 
dependent) behavior for separations > 10 and up to 
~ 1/5 of the computational box. In other words, the 
additional density fluctuations interfere with the mea- 
sured anisotropy degree at the smallest scales, making 
the statistics more isotropic. 

More interesting than confirming the scale indepen- 
dence on the structure functions in the inertial range of 
the turb ulent cascade (there is already evidence of this in 



CLV02 ) , it is to study how the degree of anisotropy de- 
pends on the turbulence parameters (sonic and Alfvenic 
Mach numbers at the injection scale). To do this we 
have calculated the degree of anisotropy on all the mod- 
els, as exemplified for model M8 in Fig. [3] and computed 
the average value in scales between 10 grid points and 
L/5. Below 10 grid points scale, the density from the 
MHD simulations is severely affected by numerical diffu- 
sion, and the effect of the large scale density gradient and 
noise are more pronounced. The results are condensed 
in Figure |4j where we have plotted the average degree of 
isotropy as a function of the sonic Mach number and of 
the Alfvenic Mach number. The Mach numbers are indi- 
cated by the different symbols (and colors in the online 
version) as displayed in the legend. The error bars show 
the variations while obtaining the average value. They 
are the product of differences in the two possible lines of 
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(a) Mean velocity (CF^x) (b) Mean velocity (CF^^) 
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Fig. 2. — Example of the iso-contours of the correlation function in one of the simulations (model M8: super-sonic and sub-Alfvenic). 
The first two panels [(a) and (b)] correspond to c orrelations in mean velocity maps, the last two [(c) and (d)] are correlations in velocity 
centroid maps (with the original density, see j ]3.1|) . In the plots on the left column [panels (a) and (c)] the LOS is parallel to the mean 13 
field. In the right column [panels (b) and (d)] mcXOS is perpendicular to the mean B field, which is in this case aligned with the horizontal 
axis. 



sight {y- and z-axis) and slight scale variations. 

One can see from Figure |4] that the centroids maps 
without noise or gradients have the smallest error bars, 
thus the least dependence on scale. The rest of the cen- 
troids and the mean LOS velocity have slightly larger 
error bars, and therefore some (small) scale dependence. 
It is also quite noticeable that, the results for the differ- 
ent centroids (with the original density or modified data) 
are very similar. This is not too surprising, because we 
have restricted ourselves to scales that were seen as less 
effected by either the cx r^^ gradient or the noise. How- 
ever, to apply this technique to real data it would be 
advisable to calculate the anisotropy degree on several 
scales and search for a scale independent range, in the 
same manner one restricts the inertial range looking for 
power-laws when analyzing power spectra. 

It is clear also, that the degree of anisotropy depends 
mostly on the Alfvenic Mach number: for increasing val- 
ues of the magnetic field the level of anisotropy increases 



as well (smaller SF±/SF\\ ). For the velocity centroids 
we have the same general trend, but not as pronounced. 
For instance, for ?J 1-4 the degree of anisotropy is 
barely distinguishable from isotropic, while in the mean 
LOS velocity map this was only the case for a very weak 
B field {Ma ^ 7.5). Centroids also show a weak depen- 
dence on Ms, but only for moderate to low magnetiza- 
tions {Ma > 0.4). For strong magnetic fields ( Ma ^ 0.4) 
we observe some dependence on the sonic Mach number. 
Since this was not noticeable in FigQa), and given the 
similarity of Figs. |4|^b-d), one can attribute such de- 
pendence to the original density field (i.e. arising from 
shocks in supersonic turbulence). This strong influence 
of the sonic Mach number is seen as a positive slope for 
the most magnetized simulations. 

The results of a clean dependence on the Alfvenic Mach 
number and a weak dependence on M^ (for small to mod- 
erate magnetic field strengths), are encouraging. While 
the sonic Mach number can be obtained by a variety 
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Fig. 3. — Example of the degree of anisotropy of the structure functions in the same model shown in Figure|2] as observed from different 
directions (a value of 1 means isotropic). The top row [panels (a)-(c)] are obtained with the mean density maps while the bottom row 
[(d)-(f)] with maps of velocity centroids. The different lines in panels (d-f) denote the density field used to obtain the centroids, the solid 
line corresponds to the original density, the dashed line to the oc gradient, and the dotted line to the addition of white noise. The LOS 
is aligned with the axis (parallel to the B field) in the left column [(a) and (d)], with the y— axis in the middle column [(b) and (e)], 
and with the z— axis in the right column [(c) and (f)]. 



of techniques ( 


Padoan et al. 1997 


Passot & Vazquez- 


Semadeni|2003 


IRowal et aim7\ 


t'ederrath et al.||2UU8 


Burkhart et al. 


2009p , the Alfvenic Mach number has re- 



along with th e Chandrasekhar-Fermi te chnique (see for 



instance Falceta-Gongalves et al. 2008) are starting to 



change this situation. 

All in all, the observable turbulent fluctuations, as rep- 
resented by velocity centroids, are sensitive to the fluid 
magnetization given by the Alfvenic Mach number Ma- 
The dependence of the anisotropy on the sonic Mach 
number Ms is, however, not always negligible. The lat- 
ter number can be ob tained using other statistics studied 
in the literature (see Kowal et al. 2007 Burkhart et al. 



2009[ |2010[ [Esquivel fc Lazarian||20lO| ) 



5. DISCUSSION AND SUMMARY 

Advances in understanding of the nature of magnetized 
turbulence drive the development of the techniques to 
study turbulence through observations. For these stud- 
ies different approaches can (and should) be used. How- 
ever, it is important to understand what is the maximal 
information that one can get from observations. 

Recently observational studies o f interstellar 
anisot ropics have been performed by [Heyer et al.| 
(|2008). The resulting anisotropics were broadly con- 
sistent with the expectations obtained in theoretical 



and numerical studies (see LPE02 EL05 ) , revealing the 
direction of the mean magnetic heMl Both theoretical 
predictions and numerical calculations, including those 
in the present paper indicate that the anisotropy ex- 
pected in the global frame of reference, which is sampled 
in the observations, should have a range in which 
they do not depend on the scale. In addition, we pay 
attention to the dependence of the velocity anisotropy 
on the sonic Mach numbers. 

In general, the measures we use to study the anisotropy 
may be affected by other interfering factors. For in- 
stance, density effects can affect velocity centroid mea- 
sures. Attempts to mitigate the effect fluctua ting density 
field on velocity centroids was attempted in Lazarian & 

modified 

However, further studie s in 
p005l); |Esquivel et al.| pOOTf 
iiled centroids" improve the re" 



Esquivel (2003) where new measures termed" 
centroids" were introduced 



Esquivel & Lazarian 
veaied that the"^^ 



modined centroids improve the repre- 
sentation of turbulent velocity only at moderate Ms . We 
have also used MVCs for the models presented in this pa- 
per and found that the results were actually worse (nois- 
ier) than with ordinary centroids. Thus in the present 
study of anisotropy we present only ordinary centroids. 

To study the robustness of the method we have 
added fluctuations to the density field and reanalyze the 
anisotropy of centroids obtained with these new data. 
We find that the range of scales in which the anisotropy 
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Fig. 4. — Degree of anisotropy in all the models averaged over scales from 10 grid points to 1/5 of the computational box. The horizontal 
axis corresponds to the sonic Mach Number, the Alfvenic Mach Number is indicated by the various symbols (and colors in the online 
version) as shown in the label. In the left panel (a) we plotted the anisotropy in velocity centroids, in the right panel (b) the anisotropy in 
the average LOS velocity. In both figures the results are obtaining averaging the two cases where the LOS is perpendicular to the mean 
field. The error bars show the maximum variation of the averaging procedure (including variation across scales). 



degree is scale independent is limited to separations 
smaller than the injection scale (in fact smaller than half 
of the injection scale), but at the same time large enough 
to avoid the effects of the additional fluctuations. In our 
models we have found that the effect of an cx gradient 
and white noise can be avoided by setting the smallest 
scale used to measure the anisotropy to 10 cells. With 
real observations one has to measure the anisotropy and 
d etermine the approp riate range. 



Heyer et al. (2008'), use not centroids, but a Principal 
Component Analysis (PCA) to study turbulence through 
observations. PCA analyzes the data in position- 



position-velocity (PPV) space which could potentially 
have some useful information that is lost in the aver- 
aging procedure in velocity centroids. However, PCA re- 
lies in calibration from numerical models, while structure 
functions and centroids can be described analytically and 
predictions of their results follow directly from theory. 

Our simulations show a dependen ce on the soni c Mach 
number which was not reported by Heyer et al. (2008j). 
We suspect the reason is that dependence on iWVis less 
prominent compared to the one on Ma- However, if 
one wants a more precise analysis it should be taken 
into account. In spite of this, we should stress the pio- 
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neering significance of the observational studies of |Heyer| 



et al. (20081 which moved this technique of magnetic 
field study fro m theoretical and numerical domain (see 
LPE02 EL05 ) to a domain of practical application. We 



believe that the technique has great future and view this 
paper as a contribution to its improvement. 

We have taken a set of simulations of fully-developed, 
driven MHD turbulence with different combinations of 
sonic and Alfvenic Mach numbers to study the velocity 
anisotropy available from observations. 

Our results can be summarized as follows: 

• Synthetic maps obtained through compressible 
MHD simulations reveal a clear anisotropy of the 
velocity field, in alignment with the direction of 
the mean magnetic field. This is true for Alfvenic 
Mach numbers less than Ma — 1.5. 

• The anisotropy, measured at scales in the inertial 
range of the turbulent cascade is scale independent 
in the global frame of reference, as opposed to the 
scale dependent anisotropy with is obtained with 
respect to the the direction of the local magnetic 
field. 



• The degree of anisotropy is dominated by the 
Alfvenic Mach number. However, a dependence on 
the sonic Mach number gets prominent for highly 
magnetized gas. The anisotropy in velocity cen- 
troids differs from that of the mean velocity. As 
many other statistics are sensitive to the sonic 
Mach number this allows to make a correction for 
Ms to retrieve AIa- 

• To show that the method is robust against large- 
scale density fluctuations and noise, we have ob- 
tained centroids with modified density fields, to 
which we have added an cx gradient, or white 
noise. We see that these have an effect that is 
mostly seen at the smallest scales. If the degree 
of anisotropy is measured avoiding such scales the 
results remained basically unchanged. 
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